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Abstract 

Treating a boundary value problem in analytical fluid dynamics, 
translation of 2D steady Navier-Stokes equations to ordinary differen- 
tial form leads to a second order equation of Riccati type. In the case 
of a compressible fluid with constant kinematic viscosity along stream- 
lines, it is possible to flnd an exact solution of the differential problem 
by rational combination of Airy's functions and their derivatives. 

1 Ordinary differential form of Navier-Stokes equa- 
tions 

Suppose to have a 2D steady flow of a fluid. Let 

$:s^(0i(s),(/>2(s)) = (x,y) (1) 

an admissible parameterization (p;) for each streamline, with $ : [a, b] 
for some suitable values a and h such that $(a) is the initial point of the 
streamline, $(6) is the end point in the considered geometrical domain. 
Then, if v = (vi,f2) is the flow velocity field, by definition of streamline 
there is a scalar function f = f{s) such that 



{VI,V2) = f{s){ij)l{s),<i)2{s)) 



(2) 



where (j) is the derivative of (j). 



Remark. The author is developing a more general model of parameteriza- 
tion for Navier-Stokes equations, in the case 3D too, under less restrictive 
hypothesis than those formulated in this work. 

Suppose that f{s) ^ 0, </>i(s) / for all s along the streamline, that is 
vi 7^ 0, and let u = v o 0. Then 

Proposition 1 Along a streamline s i — > (j){s) the following relation holds: 

V • Vvi = J-mni (3) 
Dim. Using chain rule we have 

dvi ■ dvi ■ 

But from © 

4>2 = —4)1 (5) 

Vl 



therefore 



91 

ui = — 

Vl 



dvi dvi 



(6) 

Then (jSJ follows from the fact that along the streamline ui{s) = vi{(f){s)). □ 



ox ay 



Now let cj) invertible, that is 3 g : (j){[a,b]) — > [a,b] such that s = g{(j){x,y)) 
when = Then 

Proposition 2 Along a streamline the following relation holds: 

Avi = {Vg ■ Vg)ili + iAg)ui (7) 
Dim. Prom definition, we have v = u o gi, therefore Vvi = iiiS/g. Then 
d^vi .. (dgV . (d^g\ 

— = 7/1 I I -I- 7/1 ( I 

and 



2 



therefore the thesis follows by summing the two previous expressions. □ 
The following Proposition holds for a scalar function p = p(x, y): 
Proposition 3 Let q = p o (p. Along a streamline: 

Ap = q{Ag) (10) 

Dim. From q = po (p follows p = q o g. Then, using chain rule, 

dp ^ dqdg_ ^^^^ 
dx ds dx 

and analogous equation for ^-differentiation. □ 

Now let p the density, /x the dynamic viscosity and f the body force per unit 
volume of fluid. Then the Navier-Stokes equations for the steady flow are 
(see P) 

p(v Vv) = -Vp + pf + ^Av (12) 

Identifying not derived functions /i(x, y) with their composition h o (j) and 
using simple substitutions, from previous Propositions and © we can state 
the ordinary differential form of Navier-Stokes equations: 



1 

p-puiui = pfi- q— + ^i{Vg ■Vg)ui+ p{Ag)ui (13) 
01 ox 

(t>2 

U2 = —Ui 

Note that, if the streamlines (\) = (</>!, (^2) are known, from equation ()13() it 
could be possible calculate u\. But what about the physical meaning of this 
solution? We make now some considerations about conservation of mass. 
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2 About continuity equation 



In the case of incompressible flow, the differential form of this equation is 
simply V • V = ( 1 or [Sj). In this case, from v = u o g, follows that 
dxVi = iiidxg and dyV2 = U2dyg, therefore along a streamline the continuity 
equation is 



u-Vg = 



(14) 



But it is possible, in the hypothesis (pi ^ along the streamline, to rewrite 
this equation using only ui: 

Proposition 4 //0i(s) ^ V s, the continuity equation for incompressible 
flow is 



ill + 



92 - ^(pi 

<t>l 



dg_ 

dy 



ui = 



(15) 



Dim. Differentiating the identity g((j)(s)) = s respect to variable s, we have 
Vg ■ (p = 1. Then 

9x9 = -p - —9y9 (16) 

cpi ui 

Substituting this formula in the continuity equation (flU) and using ui(p2 = 
ti2(/'i, we obtain the new identity 



ui + [4'iU2 - 4>2Ui]dyg = 
Differentiating the relation U2 = ^^^^ follows that 

01 

62^1 + (j)2Ul 4>l4>2Ul 



U2 



(17) 



(18) 



from which the relation ()15|1 is obtained eliminating 112 in ()17|1 . □ 



Remark. Suppose the streamlines are straight lines expressed by the pa- 
rameterization (j)i{s) = s, (j)2 = k, with k real constants. Then (pi = 1 and 
(p2 = 0, and from H15() the continuity equation is simply iii = 0. Then the 
velocity field is constant along a streamline, as known for incompressible 
rectilinear fiow (see e.g. jHl). 



4 



As we would investigate existence of non trivial solutions to Navier-Stokes 
equations (fT^ in the simple case x = 4>i{s) = s, we have 5x = 1 and Qy = 0, 
therefore in the incompressible case from previous Proposition follows that 
ui is constant. For a more interesting and realistic solution we assume that 
flow is steady but with spatially variable density. 

3 A steady compressible flow 

Consider a 2D steady flow where streamlines are parameterizable by the 
following expressions 

X = (j)l{s)=S, y = (j)2{s) = (j)2{x) (19) 

where (f)2 is invertible on an interval [0, L], that is for each streamline exists a 
function g : 4>{[0, L]) — > [a, L] such that s = g{x, y) = x. The function (j)2 and 
g can depend on the single streamline. Note that (pi = 1, = (1,0) and 
Ag = 0, therefore from the Navier-Stokes equations along a streamline 
become 

puiiii = pfi-q + fiili, U2 = <j)2Ui (20) 

Suppose that x-component /i of body force f is constant. Along a streamline 
we can made the realistic hypothesis that p(s) 7^ Vs; dividing the two 
members by p we obtain 

uiiii = /i - - + i/iii (21) 
P 

where = ^ is the kinematic viscosity (see ^). At this point we suppose 
that, along a single streamline, z/ and the quantity ^ are constant, with 
values depending on streamline. The latter, equivalent to p = ki J p + k2, 
can be view as a constitutive equation about the fluid. 

4 An analytical resolution 

In this section we try to find a general exact solution of the non linear 
differential equation (|21|) i . Integrating on the independent variable s and 
then dividing by i^, we have 

ni = ^ui + ^(^-fXs + ^ (22) 
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where c is an arbitrary constant. This is a first order differential equation 
of Riccati type (see or j^j). Applying the transformation 



z 



we translate previous equation into (see |2j) 



z + 



fi]s + c 



z = 



(23) 



(24) 



This is a form of Airy^s type equation (see jHI and [3j) and its general solution 
is (see |1]) 



z{s) = ciAi{t) + C2Bi{t) 



(25) 



where 



2- (26) 



(-a)3 
2zy2 



2 > . (27) 



while j4z(t) and i3i(t) are Airy's functions (see e.g. |j[|), linearly independent 
solutions of Airy's equation y — ty = which appears in optics and quantum 
mechanics phenomena. 

Using transformation ()23() , the exact solution of Navier-Stokes equation ()22|) 
along a streamline is 

m = 2u{-a) 3 — — ^ 29 

ci74z(t) + C2Bi[t) 

where Ai'{t) and Bi'{t) are the derivative of yli and Bi. 

Note that this solution has physical meaning only if a < 0. Assuming that 

inflow zone is at s = 0, usually pressure drops down towards outflow, so that 

we can assume g < 0. Also, we assume /i > 0, hence the condition a < is 

satisfied. 

Also, note that with our assumptions ui doesn't depend on the streamline 
parametric representation (x,y) = </>(s); the flow velocity field depends on 
(j) through U2 component by relation U2 = (t)2Ui- 
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5 A boundary value problem 



The general solution (PU)) depends on three constants of integration: ci, C2 
and, from l\22^ . c. Suppose we want to resolve a boundary value problem 
for equation (PT|) with ui{0) = uiq,ui{0) = uiq and ui{L) = un. The first 
two conditions, using (|22() at s = 0, give the value of c. Then, noting that 
for t G M Airy's functions Ai{t) and Bi{t) have real values (see ^7,), the 
other two constants ci and C2 can be found solving the algebraic system 

'Ul(O) = Uio, Ui{L) = UiL- 
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